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Abstract 

In this work we study by ac susceptibility measurements the evolution of the solid vortex lattice mobility under oscillating 
forces. Previous work had already shown that in YBCO single crystals, below the melting transition, a temporarily symmetric 
magnetic ac field (e.g. sinusoidal, square, triangular) can heal the vortex lattice (VL) and increase its mobility, but a temporarily 
asymmetric one (e.g. sawtooth) of the same amplitude can tear the lattice into a more pinned disordered state. In this work 
we present evidence that the mobility of the VL is reduced for large vortex displacements, in agreement with predictions of 
recent simulations. We show that with large symmetric oscillating fields both an initially ordered or an initially disordered VL 
configuration evolve towards a less mobile lattice, supporting the scenario of plastic flow. 

PACS numbers: 74.72 Bk, 74.25.Ha,74.25.Qt 

Keywords: memory effects; ac susceptibility; YBCO single crystals 



A rich variety of dc and ac dynamical behaviors aris- 
ing from the competition between pinning, thermal and 
intervortex interactions has been observed in the vortex 
lattice (VL) of type II superconductors-. Forces between 
vortices favor the formation of an ordered hexagonal lat- 
tice, in opposition to the disorder that arises from inter- 
actions with pinning centers and thermal forces, leading 
to defective polycrystalline, or glassy structures. 

In driven lattices, external driving forces participate in 
the formation of ordered and disordered structures. An 
example is the experimental 2 and theoretical 3 evidence of 
a two step depinning process of the VL as the driving 
force is increased. Initially, at low drives, the lattice un- 
dergoes plastic flow, vortices move past their neighbors 
tearing the VL and leading to the formation of a disor- 
dered lattice, with a high density of topological defects 
(e.g. dislocations). Increasing the drive above a thresh- 
old force, Ft, a dynamical crystallization occurs as pro- 
posed by Koshelev and Vinokur 3 . At larger forces vortex- 
vortex interactions dominate over interactions with pin- 
ning centers which are accounted for by an effective tem- 
perature that decreases with the VL velocity. 

Memory effects have been observed in low^£ and high 
T c materials (HTSC)^iiS, where the resistivity or the ap- 
parent critical current density J c , are found to be strongly 
dependent on the dynamical history of the VL. An in- 
crease in the mobility of the VL when set in motion by 
a temporarily symmetric (e.g. sinusoidal) ac field (or 
current is a characteristic which is common to all 
of these experiments. A proposed mechanism^ for this ef- 
fect in YBa2Cu307 (YBCO) crystals, is the annealing of 
bulk magnetic gradients in a platelet placed in a perpen- 
dicular dc magnetic field by a weak planar ac magnetic 
field. A second invoked mechanism is an equilibration 
process assisted by the ac magnetic fieldSiii. 

Transport and ac susceptibility experiments^»S^^ 
suggest that the response of a steady driven VL may 
differ qualitatively from the ac response observed in mea- 



surements involving comparable driving forces. In par- 
ticular markedly different effects for temporarily sym- 
metric or temporarily asymmetric ac drives have been 
reported 10 . It was shown that the application of a tem- 
porarily asymmetric ac magnetic field (e.g. sawtooth) 
reduces the mobility of the VL, in contrast to the effect 
of a symmetric ac field (e.g. square, sinusoidal, trian- 
gular) of similar amplitude and frequency. These effects 
were observed to be weakly frequency dependent and to 
depend strongly on the number of shaking cycles. This is 
a main result that ruled out an equilibration process as a 
possible explanation to the change in vortex mobility 10 . 
At the same time, the dynamic crystallization scenario 
becomes inadequate to describe the observed changes in 
VL mobility of the solid vortex in twinned YBCO crys- 
tals. 

Oscillatory dynamics has been described recently in 
molecular dynamics calculations'", where an oscillatory 
driving force below the crystallization threshold is able 
to order the VL after a number of cycles. The reordering 
of the VL is more efficient when vortex excursions are of 
the order of the lattice parameter ao- This VL shaking 
promotes repeated interactions between neighbors that 
heal lattice defects (i.e. the number of vortices with 5 or 
7 first neighbors). At the same time, the average mobil- 
ity of the VL increases logarithmically with the number 
of cycles. An important result is that for a tiny asym- 
metry in the amplitude of the force, the vortex lattice 
quickly disorders increasing the number of topological 
defects and the mobility is rapidly reduced after a few cy- 
cles as observed in the experiments^. On the other hand, 
when the period of a symmetrical force was increased so 
that the excursion of vortices greatly exceeded the lattice 
constant, the VL did not reorder. It was argued that if 
the period of the oscillation was large enough, the system 
should behave as when driven by steady forces below the 
threshold force. 

In this paper, we investigate the effect of shaking the 
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vortex lattice with sinusoidal magnetic ac fields starting 
from a well defined and reproducible state. In these ex- 
periments the excursion amplitude of the vortices is con- 
trolled by the amplitude of the ac magnetic field which 
is varied between 1 and 150 Oe. We found that for am- 
plitudes below a certain threshold, 20 Oe, the VL mo- 
bility increases as a function of the number of cycles of 
the shaking field. However, such an increase in the mo- 
bility is not observed above 20 Oe. Following previous 
experimental and theoretical results, we interpret this as 
an indication that plastic motion dominates the dynami- 
cal behavior of the VL. The reduced mobility would be a 
consequence of large vortex displacements which produce 
vortex lattice tearing and generation of defects. 

We measured the response of the vortex lattice to 
the shaking field by means of ac susceptibility measure- 
ments. Global ac susceptibility measurements {xac = X ' 
+ j x") were made with the usual mutual inductance 
technique in two twinned YBa2Cu307 single crystals^, 
with typical dimensions 0.6 x 0.6 x 0.02 mm 3 , and T c ~ 
92 K at zero dc field and AT C - 0.3 K (10%-90% cri- 
teria). We have obtained similar results for both crystals 
but we show the results obtained for one of them. The ac 
field, H ac , was applied parallel to the c crystallographic 
axis and the dc field, Hd C ,was applied at 8=20° away 
from the twin boundaries to avoid the effects of correlated 
defects^. In the temperature and field region of interest, 
high dissipation (x") [or low screening (x')] implies high 
mobility, and low dissipation [or high screening] implies 
low mobilityiS. 

The experiments to investigate the effect of the am- 
plitude of the shaking field followed the protocol that is 
described immediately below. First, the VL lattice was 
prepared applying the ac configuration field (H c /) for a 
number of cycles (N ~10 5 , / = 10 kHz). A temporarily 
asymmetric (sawtooth) H c / was used to prepare a low 
mobility configuration, LMC. A temporarily symmetric 
(sinusoidal) H c / was used to prepare an initial vortex 
configuration with high mobility, HMO It is worth to 
note that the lattice was always prepared in the LMC 
before applying the symmetric H c / to attain the HMO 
The configuration field and the measuring H ac field were 
provided by the susceptometer primary coil. After set- 
ting the desired starting configuration, H c j was turned 
off and a temporarily symmetric shaking ac field (R s h) 
with an amplitude that could be varied between 1 Oe 
and 150 Oe (/ = 1 Hz) was applied. This field was sup- 
plied by the same electromagnet that provided the dc 
field. After shaking the vortex lattice for a number of 
cycles (N s /j), H s h was turned off and the first harmonic 
of the magnetic ac susceptibility was measured with a 
sinusoidal ac field with an amplitude H ac = 2 Oe and / 
= 10.22 kHz. The measured susceptibility is related to 
the mobility of the VL and quantifies the effect of the 
shaking field in it. With this protocol we can study the 
effect of ~N s h cycles of the shaking field to an initial vor- 
tex configuration with either low (LMC) or high (HMC) 
mobility. This procedure is repeated for each measured 



value of N s /,. 

Figure 1 shows ac susceptibility measurements x" vs 
T cooled in dc (H dc = 2 kOe, O = 20°) and ac fields 
(H QC = 2 Oe, / = 10.22 kHz ). Differences between VLs 
with high and low mobility are measured for tempera- 
tures below the melting line (solid VL). We choose to 
make our measurements at T^85 K where a larger dif- 
ference between the measured signal of the high and low 
mobility states is observed. The lower point (point A) 



FIG. 1: Ac susceptibility measurements x" vs T. This mea- 
surement was made lowering temperature with H ac and Hd c 
turned on. 



was obtained after field cooling to 85.5 K, turning off the 
measuring field and applying 10 5 cycles of a sawtooth 
ac magnetic field, H c / = 7.5 Oe and / = 10 kHz. The 
7.5 Oe ac field (that penetrates the sample completely^) 
was turned off and ac susceptibility was measured. The 
application of the temporarily asymmetric drive reduced 
vortex mobility, and we call this a low mobility configu- 
ration, LMC. The higher point (point B) was the dissi- 
pation level obtained after setting a LMC (point A) and 
then applying 10 5 cycles of a sinusoidal ac field, H c / = 
7.5 Oe and / = 10 kHz. The 7.5 Oe ac field, was turned 
off and the ac susceptibility was measured. The mobility 
of the VL is clearly enhanced as a result of the applica- 
tion of the sinusoidal field. We call this vortex state a 
high mobility configuration, HMC. 

Our next experiments were performed at fixed temper- 
ature and dc magnetic field. As anticipated above, the 
temperature was chosen to correspond with the low tem- 
perature maximum in V' (see Fig. 1) in order to obtain 
a large observable difference between the measured sus- 
ceptibility for low and high mobility VL configurations. 
For our samples, for a dc magnetic field of 2 kOe and 
with our selected measuring ac field, this temperature is 
around T= 85.0 K. 

In Fig. 2 we show x" vs the number of cycles (N s h) of 
the shaking field. The shaking field was chosen to have 
different amplitudes (from 8 to a 120 Oe) and was applied 
to a LMC prepared with H c / =7.5 Oe. Equivalent results 
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FIG. 2: x" vs number of cycles of the shaking field (N s h) for 
different amplitudes, starting from a low mobility configura- 
tion., LMC. 



FIG. 3: x" vs number of cycles of the shaking field (N 3 /j) for 
different amplitudes, starting from a high mobility configura- 
tion HMC. 



are obtained for x'- For shaking fields with the lowest am- 
plitude (8 Oe), x" increases roughly as the logarithm of 
N s /,. The same dependence is observed for intermediate 
amplitudes (10-20 Oe) for which it is also observed that 
x" increases with the amplitude of the shaking field (for a 
given number of cycles, N a /,). It is interesting to note that 
this increasing trend in x" does not continue if the am- 
plitude of the shaking field is further increased. Instead, 
we observed that x" reaches a maximum at around 20 
Oe and starts decreasing for larger amplitudes. For H s h 
above 80 Oe, the x" of the final state (after 1000 cycles) 
is in fact comparable to the x" of the initial LMC. 

As a larger x" implies a larger VL mobility, these re- 
sults indicate that there is an "optimum" amplitude of 
H s /j (~20 Oe) for which a maximum mobility in the VL is 
obtained (for a given number of cycles). They also show 
that a high amplitude symmetric ac field is not effective 
in reordering and increasing the mobility of an initially 
disordered VL. By comparison with numerical simula- 
tions, it appears that large vortex displacements produce 
plastic tearing of the VL and the overall response be- 
comes equivalent to the response of a lattice in a low 
mobility configuration. 

We also studied the effect of the shaking field on an 
initially ordered VL (HMC). In Fig. 3 we show x" vs 
the number of cycles of the shaking field (N s h) of differ- 
ent amplitudes starting from the HMC. At low amplitude 
magnetic fields (<10 Oe) x" (and the VL mobility) stays 
approximately constant but at higher amplitudes (20-80 
Oe) an overall reduction in x" is observed. Even one cy- 
cle is enough to significantly alter the VL configuration. 
As N s h is increased, the value of x" seems to go through a 
shallow minimum but the x" never recovers to the initial 
value at N s ^=0. Shaking fields with an amplitude larger 
that 80 Oe strongly reduce x" to values which are com- 
parable to the one in an LMC. As discussed above, the 
results in Fig. 2 show that a large amplitude symmetric 



magnetic field is not effective in reordering and increas- 
ing the mobility of the VL. Fig. 3 shows that such an 
oscillating magnetic field also distorts an initially ordered 
VL and reduces its mobility. 

The effects of the amplitude of the shaking field H s /j 
are more clearly observed in Fig. 4 which shows x" as 
a function of the amplitude of the shaking field for a 
fixed number of cycles(N s / l =200). We show measure- 
ments that were performed starting with high and low 
mobility configurations. Starting with a LMC it can be 
seen that the ability of the shaking field to improve VL 
mobility increases up to a maximum (K s h ~ 10-20 Oe) 
and then decreases. For amplitudes H s /j ^80 Oe the 
VL seems to have reached a configuration just slightly 
different from the starting LMC, i.e. a high amplitude 
sinusoidal field does not remove VL defects. For a HMC 



FIG. 4: x" vs amplitude of shaking field for N s h=200, from 
Low and High mobility configurations. 



starting state, low amplitudes do not modify the dynam- 
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ics of the VL. As H s /j is increased above 20 Oe, there 
is a clear reduction in mobility, and higher shaking am- 
plitudes configure the VL close to the LMC state. In 
fact for H s h > 40 Oe the final mobility is independent of 
the starting configuration. It is interesting to note that 
this result is independent of the frequency of the shaking 
field in the range tested (0.1 Hz < f s h < 3 Hz). Given 
that the dissipation in the sample is directly related to 
the number of cycles per unit of time, this result implies 
that the observed effects are not related to local heating. 

Following Ref.7, when the average vortex excursion 
produced by the shaking field is comparable to the lat- 
tice constant, ao, the VL mobility increases as the vortex 
lattice orders and moves in an increasingly coherent way 
(the calculated number of defects in the lattice and its 
mobility vary as the logarithm of the number of cycles of 
the oscillating force). On the contrary, calculations pre- 
dict that when the excursion of vortices greatly exceeds 
the lattice constant, the VL mobility is reduced as the 
plastic motion tend to increase disorder. In order to re- 
late our results with theoretical predictions, we estimated 
the average vortex displacement under the action of the 
oscillating ac field. The distance that a vortex at the 
sample boundary moves when a field perturbation H s /j is 
applied can be roughly estimated by: 



where r is the sample radius, assuming the magnetic in- 
duction B^H^c- For our experimental conditions, r^0.3 
mm and Hdc — 2000 Oe, and considering a triangular 
lattice (ao ~0.1/xm), <u> ^ao occurs at a shaking field 
amplitude H s h ~2 Oe. However, in our experiments the 



maximum dissipation (implying maximum mobility) oc- 
curs for H s /j ~10 Oe (see Fig. 4), so that the above 
approximations slightly underestimate the shaking field 
limit. Note that the exact value at which the maximum 
in x" is observed could depend on the rigidity of the vor- 
tex lattice and the density of pinning centers. The more 
rigid the lattice the more difficult is to create defects in 
it. This implies that the field amplitude estimated above 
for the position of the peak is a lower limit and could 
increase with vortex rigidity. 

We find that our results are qualitatively in accordance 
with numerical simulations indicating that if vortices are 
forced to oscillate with amplitudes larger than the typ- 
ical VL parameter, plastic motion introduces topolog- 
ical defects and reduces mobility. For smaller drives, 
vortices perform small excursions interacting repeatedly 
with neighbors and as the lattice becomes successively 
more ordered, its mobility increases. 

To conclude, in this paper we have shown that tem- 
porarily symmetric vortex oscillations, forced by sinu- 
soidal ac fields, increase the mobility of the VL in twinned 
YBCO crystals. However when the amplitude is larger 
that a certain threshold, the temporarily symmetric os- 
cillation reduces the mobility. This is an indication that 
large vortex displacements may produce vortex lattice 
tearing and the mobility is reduced. We have also shown 
that a healed lattice with initial high mobility can be 
torn by the driven symmetric oscillation, if vortex dis- 
placements are much larger than the VL parameter. The 
ordered VL reduces its mobility even with just one oscil- 
lation of the shaking field. 
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